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We study the kinetic equilibration of gluons produced in 
the very early stages of a high energy heavy ion collision in a 
"self-consistent" relaxation time approximation. We compare 
two scenarios describing the initial state of the gluon sys- 
tem, namely the saturation and the minijet scenarios, both 
at RHIC and LHC energies. We argue that, in order to char- 
acterize kinetic equilibration, it is relevant to test the isotropy 
of various observables. As a consequence, we find in particu- 
lar that in both scenarios elastic processes are not sufficient 
for the system to reach kinetic equilibrium at RHIC ener- 
gies. More generally, we show that, contrary to what is of- 
ten assumed in the literature, elastic collisions alone are not 
sufficient to rapidly achieve kinetic equilibration. Because of 
longitudinal expansion at early times, the actual equilibration 
time is at least of the order of a few fermis. 

PACS numbers: 25.75.-q, 12.38.Mh 



I. INTRODUCTION 

The main interest of studying very high energy heavy- 
ion collisions is the possibility that a quark-gluon plasma 
(QGP) be formed. It is widely believed that a large 
amount of real partons (essentially gluons) with trans- 
verse momentum p t ~ 1 — 2 GeV is freed in the very early 
stages of such collisions, carrying most of the produced 
transverse energy Jl|,||. While rapidly diluted by longi- 
tudinal expansion, this parton gas tends to locally equi- 
librate via mutual interactions between its constituents. 
When the mean free path of a parton becomes too large, 
the description in terms of partonic degrees of freedom 
ceases to be valid; one is instead left with an expanding 
gas of hadrons. Whether the dense parton gas reaches a 
state of kinetic equilibrium before it hadronizes is a ques- 
tion of great interest for the interpretation of the data at 
RHIC and LHC. So far, most of the calculations con- 
cerning experimental signatures of the QGP rely on the 
assumption that the system is at least kinetically equili- 
brated. 

It was commonly assumed for a long time (see e.g. 
that shortly after the nuclei collide, kinetic equilibrium 
is rapidly achieved on time scales < 1 fm (see e.g. ||). It 
was however shown that, due to the effect of expansion, at 



least elastic collisions may not be effective enough 
An important question is also to characterize the initial 
state of the parton system formed just after the colli- 
sion ||,|7j . The purpose of this paper is to reexamine the 
role of elastic collisions in the process of kinetic equilibra- 
tion, assuming that, already at early times, a local Boltz- 
mann equation may be written for the partonic phase 
space distribution. We specifically compare within a sim- 
ple relaxation time approximation the case where initial 
conditions are given by the saturation scenario, which 
may already be valid at RHIC, to the case where the ini- 
tial gluons are produced incoherently in the perturbative 
QCD regime, the so-called minijet scenario. We consider 
only the pure gluon case and 2 — > 2 elastic collisions in 
the limit of small- angle scatterings. By calculating the 
time-dependent relaxation time in a self consistent way, 
we are able to obtain various observables as functions of 
time, thus probing kinetic equilibration. The method is 
described in sections [n] to In section we present 
our results. Comparing with Ref. where the exact so- 
lution is worked out in the saturation scenario, one can 
assess the validity of our approach, which appears to be 
surprisingly good. We then apply this method to the 
minijet scenario, considered in Ref. |}|, and we consis- 
tently take care of conservation laws. We argue that, to 
characterize equilibration, a reliable criterium is to test 
the isotropy of various observables. As a consequence we 
find in particular that for both initial conditions the sys- 
tem does not reach kinetic equilibrium at RHIC energies. 
We discuss the reliability of our description and argue in 
particular that, due to the fragility of the weak coupling 
approximation, it appears difficult to obtain definite con- 
clusions at LHC energies. In any case, the kinetic equi- 
libration time is an order of magnitude bigger than the 
typical 1 fm estimate usually assumed. The hypothe- 
sis that elastic collisions are sufficient to rapidly achieve 
kinetic equilibration, often made in the literature j3| is 
thus shown to be wrong, as already suggested in [0 . Sim- 
ilar conclusions have been obtained in a slightly different 
context by the authors of Ref. |10l . 
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II. THE BOLTZMANN EQUATION IN THE 
SMALL ANGLE APPROXIMATION 



Shortly after the gluons have been produced, the oc- 
cupation number becomes small enough so that one can 
treat the further evolution by means of a Boltzmann 
equation describing the time evolution of the local dis- 
tribution f(p, x, t) for on-shell particles^ When the oc- 
cupation number becomes smaller than unity, one can 
make the further approximation of a gas of classical par- 
ticles. Here, we consider only the gg — ► gg elastic pro- 
cesses in the collision term, which we shall treat in the 
leading small-angle scattering approximation. Assuming 
one-dimensional expansion in the early stages after the 
collision and longitudinal boost invariancc in the central 
region in the center of mass frame, the Boltzmann equa- 
tion can be written in terms of a derivative at constant 
P.t fill: 



d t f(f,t)\^ t =C(p,t). 



(2.1) 



Defining the moments N s = (p s ) — J- p s f(p, t), with the 
notation J~ = 2(N 2 — 1) J d 3 p/(2n) 3 , one obtains for the 
collision term C, in the limit of small- angle clastic scatter- 
ings between classical particles (/ <C 1), to logarithmic 
accuracy p|,|T3[|, 

C(p,t)=LN V 2 p f(p,t) + 2LN- 1 V p [tT •/(?,*)] , (2.2) 



p/p, and where C 
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where v 

Xmin being the minimum scattering angle, whose value 
is determined by the physics of Debye screening. We shall 
come back to this later. 

In the present paper, we will solve the Boltzmann equa- 
tion in a relaxation time approximation, that is we will 
solve the equations for moments of the distribution /, 
taking into account howe ver the microscopic dynamical 
information contained in (2.2). 



III. THE "SELF-CONSISTENT" RELAXATION 
TIME APPROXIMATION 



-(/ - f eq )/e, with f eq (p,t) = A(t) exp(-p/T(t)). We 
assume that the relaxation time 8 does not depend on 
p. This has to be viewed as a mean-field like approx- 
imation : the different modes are effectively decoupled 
and their evolution is governed by an effective relaxation 
time which contains the microscopic information and is 
to be computed self-consistently. Indeed, strictly speak- 
ing this equation has no solution and one has to take it 
in a weaker sense, that is at the level of moments of the 
partonic distribution: 



m(p,t)C(p,t) = - 



(m)(t)-(m) eq (t) 

e m (t) 



(3.1) 



where m(p,t) is an arbitrary function, (m)( eg ) = 
Jp m(fi, t) /( eg ) (p, t), and 6 m {t) is the associated relax- 
ation time, different choices for the function m leading 
to different relaxation times. 

Let us be more precise on this last point: it is clear that 
the RTA, being a one-time-scale ansatz, cannot fully re- 
produce the evolution of the whole ensemble of modes. 
A particular choice for the function m picks up a partic- 
ular momentum scale in the momentum distribution /, 
whose contribution to the moment (m) dominates. The 
relaxation toward isotropy of this part of the distribution 
is characterized by the time scale m . So the appropri- 
ate choice for the moment function m depends on the 
momentum scale characterizing the physics one wants to 
study. 

In this paper we are concerned with the question of 
kinetic equilibration, that is with the isotropy of "ther- 
modynamic" quantities, like longitudinal and transverse 
pressures Pt,(t) = (p 2 /p) and Prit) = (p]_/p), where 
Pj_ = (Px +Py)/2- The particles we are interested in have 
energy of the order of the average energy per particle 
e(t) = e(t)/n(t), (n — Nq and e — N\ are the parti- 
cle number and energy densities per unit volume respec- 
tively). A pertinent choice leading to simple equations 
is^| m(p) — (p 2 z —p±)/p or equivalently, (m) = Pl — Pt- 
Inserting the Landau-Mueller collision term fl2.2| ) in (3.1) 
and introducing the notation^: iV* = (p'j, p s ~ 2 ) and 
N^- = (p 2 L p s ^ 2 ) 1 one obtains the following equation for 
the relaxation time: 



The Relaxation Time Approximation (RTA) simply 
consists in replacing the effect of the collisions in the 
Boltzmann equation (2.1) by an exponential- like relax- 
ation towards a local "equilibrium" distribution^]: C = 



N z — N 1 - 

1 Q 1 = 4£ N (JV*j - N±J + 2 C JV_i (7V Z - Nq) . 

(3.2) 

The two other parameters of the relaxation time 
ansatz, X(t) and T(t), are computed using conservation 



1 One can argue that the requirement for using a local Boltz- 
mann equation for on-shell quanta is that the occupation 
number / be < 1/as, where as is the strong coupling con- 
stant |ll|. 

2 The distribution f eq has the meaning of a local equilibrium 
distribution only when the system has reached the hydrody- 
namic regime, where A(t) and t 1 ^ 3 T(t) are constants |14j 



The choices (m) = Pl or (m) = Pt give the same results 
as those presented below. The choice (m) — e leads to the 



equation of energy conservation (see Eq. (3.3)) which contains 
no information about collisions. 
4 For example, P L = N?, Pt = . 
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laws: conservation of energy, and of particle number in 
elastic collisions^]. In the RTA, these give the relations 



.Ni - 1 



<*) = e eq (t) = 6 y2 
n{t) = n eq (t) = 2 ' ""-r-' 



X(t)T 4 (t), (3.3) 
., X(t)T 3 (t). (3.4) 



Note that the particle number conservation law implies 
that the number density exactly falls like tn(t) = 
to n(to), with to the time when our description begins. 



IV. SOLVING THE EQUATIONS 

The solution of the Boltzmann equation in the RTA is 
given by p) 

f(puP z ,t) = f (pt,p z ^)e- x U 
to 

t x(t')-x(t) + 

dt' - 



9{t 



— fe q (pt,Pzj;,t'), (4.1) 



a se cond estimation of M(ti + St), again using 

Eq. (4.2). We repeat these steps until convergence of 
our three parameters to some fixed accuracy. We have 



checked the validity of this method by solving Eq. (4.2) 
in some exactly solvable cases, for example by com puting 
n = No at each time-step with the help of Eq. (4.2), with 



^ 0) (a) 



V. INITIAL CONDITIONS AND RESULTS 

In this section, we apply the above formalism to the 
study of kinetic equilibration of the gluon gas formed 
during the very early stages of the collision. We study 
two different scenarios for the initial condition, namely 
the saturation and minijet scenarios. The corresponding 
initial distributions are taken from Ref. Q for the former 
case (see also j|) and from Ref. || for the later. 

To proceed further, one has to give a prescription for 
regularizing the logarithmic singularity C appearing in 
the collision term (Eq. (2.2)), in the small angle approx- 
imation. 



where x(t) — J^dt' 8 1 (t'), fo(p) is the initial distribu- 
tion, and where X(t), T{t) and are computed at 
each tim e with the help of Eqs. ( |3.2| )-(3.4). We can use 
Eq. (4T) to give an integral expression for all the quan- 
tities appearing in these equations. In the present case, 
every phase space integral can be computed analytically, 
thus considerably simplifying the numerical task. One 
typically obtains a formula of the type 



M(t) = M(to)T^(to/t)< 



+ j t dt' e —T^\t'/t)M eq {t'), (4.2) 

where x = x(t), x' = x(t'), 9' = 6{t') and M = N S , N% or 
(s = — 1, 0, 1), and where the corresponding functions 

The function depends 



Fj^ an d ^m" are calculable. 



r(e<?) 



on the form of the initial distribution. The "equilibrium" 
moments M eq , computed with f eq , only depend on A and 



T (see Eqs. (Al) and (A2) in the appendix} 



We solve the system of equations (p.2j)-(|3.4|) and (|4.: 
by using a method similar to that described in Ref. 
suppose we know the parameters A, T and Q" 1 for ev- 
ery time t < t±. We want to compute their values at 
time ti + St. We use a first guess (for example the 
values at time ti) to compute a first e stim ation 
of the moments M(ti + St), using Eq. (4J2). We then 
compute a second esti mat i on o f our three parameters 
with the help of Eqs. (3.2)-(3.4), with which we obtain 



5 The "fugacity" parameter X(t) is needed to enforce particle 
number conservation in elastic collisions. 



A. Screening 

The logarithmically divergent integral appearing in the 
small scattering angle collision integral is physicall y regu - 
lated by screening effects in the gluon medium |^, |l3| , |l5| . 
The minimum scattering angle is given by the relation 
m 2 D = p 2 Xmim where mu is the Debye mass, an expres- 
sion of which can be obtained in terms of the distribu- 
tion function in the linear response approximation p5j , 
and where p is the typical momentum of particles in the 
medium. We shall see below that for early times, because 
of the expansion, the particles in the central region have 
essentially zero longitudinal momentum]^]. In this situ- 
ation, the exchanged gluon is essentially transverse and 
the relevant screening mass is the transverse mass B 



i 2 T (t) = 



a s N c 



— fM = l^l N - l{th 



and we write L 



2 ■ ) 

nun J 



Pt 



Ps 



P 



2 and where 



ln((p 2 ) / m 2 ") , where 
.) /n denotes the av- 
erage per particle. Although the above choice for L is 
motivated by the highly anisotropic form of the distri- 
bution at early times, we shall use it also when the sys- 



This is explicit in the initial dis trib ution corresponding to 
the saturation scenario (see Eq. (5.1)). In the minijet sce- 
nario, although the initial distribution is isotropic in momen- 
tum space, it becomes rapidly peaked around p 2 = (see 
Fig. |). 



3 



tern approaches kinetic equilibrium^. It should be noted 
here that the validity of the logarithmic approximation is 
marginal because of the fact that the coupling constant 
is not very small compared to one (in particular, one has 
typically asN c ~ 1, see below)|^. 



B. Saturation scenario 

In the saturation scenario one assumes that the pro- 
duced transverse energy per unit rapidity is dominated by 
small-x gluons having a transverse momentum p t ~ Q s , 
the saturation momentum ||. These gluons are essen- 
tially freed during the collision, after a time U ~ 1/Q S 
which is also the time needed for gluons in different units 
of rapidity to physically separate from each other. Af- 
ter a time U, gluons populating the central region in the 
center of mass frame of the collision have essentially no 
longitudinal momentum (p z ti ~ p z /Qs ^ !)■ Following 
Mueller §], one writes for the occupation number in the 
central region at U: 



fsat(P,U) 



a s N t 



5{PzU)Q{Q 2 s 



where c ~ 1 is a numerical constant, which has been com- 
puted numerically using a classical field approximation in 
Rcf. @. For an SU{2) gauge theory, c = 1.3. The ap- 
proximation of classical particles makes sense when the 
occupation number / < 1 so that, following the authors 



of Ref 



we take, as initial time, t 



(i 



a s AT, 



U where 



one has, assuming free streaming between ti and to, 

Hp) = fsat(p,t ) = 6(pA) &(Q 2 S -Pt) ■ (5.1) 

The values of the parameters corresponding to RHIC and 
LHC energies given in || are summarized in Tab. |[ In 
the following we take as = 0.3, N c = 3. We evolve 
the system until the particle number density becomes 
lower than 1 per fm 3 , when the description in terms of 
partonic degrees of freedom becomes meaningless. The 
corresponding time t max is evaluated using the conserva- 
tion of the particle number: t n(t) — t n(i ) (see Tab. |). 
This rough estimate serves only to give an upper limit for 
our description. In particular one obtains t max ~ 30 fm 
at LHC, which is too long for the approximation of lon- 
gitudinally boost invariant geometry to be valid. Never- 
theless, we present our results for t < t max in order to 
show the behavior of the system. 

Fig. Q shows the logarithm L as a function of time. 
At early times, L <C 1 whatever the way one chooses 



7 For an isotropic distribution, a more appropriate choice 
would be Li ao — ln((p 2 ) / 'm 2 D ) . In this case however m 2 D = 
2ml, (Pt) = 2 (p 2 )/3, and L - L iao = In 4/3 ~ 0.3. 

8 For more details on this point, see |IL6[. 



to implement screening. As time goes on, although the 
situation becomes better, the logarithmic approximation 
remains marginal. 

We first test our relaxation time approximation by 
comparing our results with those of Ref. We show 
in Fig. H the time evolution of the energy density e and 
the longitudinal and transverse pressures Pl and Pt at 
RHIC, as well as those of the mean longitudinal and 
transverse squared momentum per particle (pi) and (p 2 ^), 
both for RHIC and LHC energies. These are to be com- 
pared with Figs. 7 and 9 of ||] (we choose the units so 
as to make the comparison easiest): one observes a semi- 
quantitative agreement with the exact solution, which 
is remarkably good in view of the simplicity of our ap- 
proach. 

We now come to the study of kinetic equilibration of 
the gluon gas. To this end, we measure the anisotropy of 
the distribution by means of the ratios 



Rk 



(Pl/P k ) 
{P'i/P k 



N z 



(5.2) 



which should approach 1 as the system approaches ki- 
netic equilibrium. We stress here that the criterium we 
use to characterize kinetic equilibrium (the isotropy of 
the distribution) is quite different from that used in ||. 
As a consequence, although we have essentially similar 
results, we arrive at quite different conclusions. Indeed, 
the time evolution of Rk for k — 0, . . . , 3, at RHIC 
and LHC energies, is shown in Fig. ^, where one can 
see that the system is far from being isotropic even at 
t ~ 10 fm at RHIC. The situation looks better at LHC, 
where Rk — 0.8 for t ~ 10 fm. It is to be noted that in 
this scenario, the only scale is Q s and consequently the 
results for RHIC and LHC are essentially the same when 
everything is expressed in units of Q s . 



C. Minijet scenario 

In the minijet scenario, the initial gluons are produced 
by means of hard and semi-hard collisions between the 
partons of the incident nuclei, and one can compute the 
initial energy and particle number density in perturbative 
QCD jij. The initial gluon multiplicity is dominated by 
the softest of these hard gluons, having an energy po ~ 
1 — 2 GeV, which are produced on a typical time to ~ 
1 /pq. Following the author of Ref. |j| , we parametrize the 
initial state by a Boltzmann distribution. In the central 
region, we have 



fo(f) = fjet(p,t = 1/po) 



-P/Tj. 



(5.3) 



where the parameters Xj e t = \{to) and Tj et — T(to) are 
determined from the initial energy and particle number 
densities: e jet = e{t ) = e eq (t Q ) and n jet = n(t ) = 
n eq (to). We recall the values of the parameters used in U 
for RHIC and LHC in Tab. 0. As in the previous sec- 
tion, we use as — 0.3, N c = 3. Concerning the validity 
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of the logarithmic approximation, the situation is more 
confortable here: in the time interval we considered, we 
have 1.5 < L < 5. 

Fig. [| shows the time evolution of the ratio Rk de- 
fined in ( jo~2] ) for both RHIC and LHC energies. As ex- 
pected, one observes a regime a free streaming at very 
early times, until the effect of collisions begins to set in, 
that is for time t — t <C 9(t ). Concerning the degree of 
isotropy the system reaches, the situation is even worth 
than in the saturation scenario: the system is always far 
from isotropy at RHIC; at LHC, the system is still quite 
anisotropic (Rk ^ 0.8) until t ~ 30 fm. Here again, we 
note that the only scale of the description is Tj et , and 
the results for RHIC and LHC are essentially the same 
when expressed in units of Tj e t . 

The above conclusions are to be compared with those 
of Ref. ||, where the authors obtain relatively short 
"equilibration times" , of the order of 4 — 5 fm, both at 
RHIC and LHC. However, although only elastic collisions 
between gluons are considered, the total number of par- 
ticles is not conserved in Q , and in fact grows with time 
(the number density n falls slower than 1 /t) . This means 
in particular that collisions between particles are anoma- 
lously frequent and consequently the system equilibrates 
anomalously fast. The similarity of equilibration times 
obtained in || at RHIC and LHC is due to the fact that 
the authors work with a running coupling constant which 
grows with time, as the typical momentum of the parti- 
cles decreases. The case of a running coupling constant is 
examined in next section. We shall see that our previous 
conclusions remain unchanged. 

D. Discussion 

Let us first remark that the model used here to describe 
the evolution of the gluon gas can only give qualitative 
information. Here we use it to study the question: are 
the elastic collisions efficient enough to achieve kinetic 
equilibrium, and if so on which time scale? In order to 
obtain a reliable answer, let us study the sensitivity of 
the above results to the details of the description, for 
example the value of as or the prescription used to cut- 
off the logarithmic divergence in the small- angle limit. 

1. Running coupling constant 

It is interesting to consider the case where the coupling 
constant is allowed to run, growing with time as the typ- 
ical energy of the particles in the medium decreases. As 
mentionned above, this should accelerate the equilibra- 
tion process. Here we shall compute the coupling con- 
stant at each time as a function of the mean energy per 
particle: as{t) = ag(/x = e = 3T), with 

asM=as(M z ) , (5-4) 



where Mz — 90 GeV the mass of the Z° boson, 
a s {M z ) = 0.1, and where we take A QCD ~ 200 MeV. 

In the minijet scenario, at LHC, we have 0.2 < ag(t) < 
0.4 on the time interval we considered, and we obtain very 
similar results to the case where as = 0.3. At RHIC, 
0.3 < as(t) < 0.5 and, although the situation is better 
when compared to the fixed as case, the system still does 
not reach kinetic equilibrium. The corresponding curves 
for the ratio Ri = Pl/ Pt are shown in Fig. [|. 

In the saturation scenario, the mean energy per parti- 
cle is smaller than in the previous case and the coupling 
constant becomes rapidly > 0.5. and the results have no 
physical significance. This is to be viewed as a manifesta- 
tion of the fragility of the weak coupling approximation. 

2. Robustness of the results 

In order to estimate roughly the accuracy of our re- 
sults, we give two estimates, simply replacing L — > 2L 
and L — ► L/2 in the collision term, keeping as = 0.3. 
The corresponding curves for the ratio R\ — Pl/ Pt are 
shown in Figs. [| for both scenarios, at RHIC and LHC 
energies. In the minijet scenario we also show the curves 
corresponding to the case where the coupling constant 
depends on time, which are seen to be contained within 
our "uncertainty bands" . One observes that the quali- 
tative conclusions of previous sections are unchanged: in 
both scenarios, elastic collisions are not sufficient to drive 
the system towards kinetic equilibrium at RHIC energies. 
At LHC, kinetic equilibrium may be reached in the sat- 
uration scenario even for times reasonnably short for the 
approximation of longitudinal expansion to have a sense, 
say t ~ 10 fm. However, it appears difficult to draw 
any firm conclusion because of the uncertainties of our 
qualitative description. In the minijet scenario at LHC 
energy the situation is even worst and one cannot draw 
definite conclusions concerning kinetic equilibration even 
for very long times. What one can say however is that, in 
both scenarios, the assumption that kinetic equilibration 
is rapidly achieved by means of elastic collisions only, on 
time scales < 1 fm, is not reliable at RHIC energy and 
questionnable at LHC energy, where the actual equilibra- 
tion time is at least of the order of a few fcrmis. 



VI. CONCLUSION 

We have studied the kinetic equilibration of the gluon 
gas initially produced in very high energy nuclear colli- 
sions, by considering 2^2 small angle elastic scatter- 
ings only. Two different scenarios have been assumed for 
the production of such gluons: the saturation and the 
minijet scenarios. By using a simple "self-consistent" re- 
laxation time approximation, we are able to reproduce 
semi-quantitative features of the exact solution of the 
Boltzmann equation. By measuring the anisotropy of 
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different observables as a function of time, we can follow 
the system toward kinetic equilibrium. Our results show 
that elastic collisions are not as efficient as usually be- 
lieved to achieve kinetic equilibrium. This contradicts a 
widely used assumption. 

This does not means that kinetic equilibrium is not 
achieved in a real collision. Indeed, concerning the phys- 
ical question of kinetic equilibration, the present study 
is only preliminary as one should include other contri- 
butions in the collision integral. In particular, inelastic 
branching processes have recently been proposed as play- 
ing a very important role in this respect M. 
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Eqs. (|Al| )-( A2) can be written in the generic form of 
Eq. (|4.2| ), using the following definitions: T e ^ (a) = 
ah s (a) and ^"^(a) = 3a 3 h z s (a). 
Writing A = \/\ — a 2 , one obtains: 



arcsinA a + /i_i(a) 
1 , = 



arctan A 2 
h-2(a) = -. , h 2 {a) = 



2 

,-,2 



A 3 

The functions h z s can be deduced from the previous ones 
by using the obvious relation 

For the particular case s = 0, one obtains 

1 / arctan A \ 
A 2 ^ A ) ' 



h (a) 



1 



K(a) 



APPENDIX: MOMENTS OF THE 
DISTRIBUTION 



We give here the explicit form of the equilibrium mo- 
ments M eq and of the functions T^/^ appearing in 



Eq. ( |4.2| ), for the different moments M used in this 
paper. These are of three different sort: "isotropic", 
as N s — (p s ), "longitudinal", as N* — {p 2 p s ~ 2 ), and 
"transverse", Nf_ = (p\p s ~ 2 ), with the obvious relation: 
N s = N* + 27V 1 - . Introducing the notation 

XAWm - / A if = feq(jP,to) 

{A\m-< B i£f o (0)=6(p,)g(p ± ), 



one easily gets, using Eq. (|4.1[) , 

te x N s (t) =t N s (t ) {h s (t /t)\\l} 
f* dt' 

+ 1 -Qft'e* Nt"{t')h s {t'/t), 

and 

t 3 e x N:(t)=t 3 Q N s (t ) {K(t o /t)\\0} 



(Al) 



where 



+ / ^(t'f^ N%«(t')K(t'/t), (A2) 



N 2 -l 



Nt«(t) = p s f eq (p, t) = (s + 2)\ A(t) T s+i (t) , 



and where 



h s (a) = / dx [1 - (1 - a 2 )x 2 ] ' , 
Jo 

„2 n ft „2n„21«/2-1 



h'Ja) 
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SATURATION 


Qs (GeV) 


to (fm) 


n(to) (fm" 3 ) 


e(to) (GeV/fm 3 ) 


tma: (fm) 


RHIC 


1. 


0.4 


18.1 


12.0 


~ 10 


LHC 


2. 


0.18 


163.4 


217.9 


~ 30 



TABLE I. Values of the parameters Q a and to characterizing the initial state in the saturation 
scenario at RHIC and LHC energies. The energy and particle number densities are computed from 
the initial distribution Eq. (^l]) . The time t max is the time after which the particle density becomes 
less than 1/fm 3 and is evaluated using particle number conservation: tn(t) — ton(to). 



MINUET 


to (fm) 


rijet (fm 3 ) 


e jet (GeV/fm 3 ) 


A jet 


T jet (GeV) 


tmax (fm) 


RHIC 


0.18 


34.3 


56.0 


1.0 


0.535 


- 10 


LHC 


0.09 


321.6 


1110.0 


1.0 


1.13 


~ 30 



TABLE II. Values of the energy and particle number densities characterizing the state of the 
gluon system at the initial time to in the minijet scenario, at RHIC and LHC ener gies. The values of 
the parameters A le t and Tj et characterizing the assumed initial distribution Eq. fl5.3|) are obtained 
from Eqs. ([j.3])-(3.4) at t = to- The time t max is the time after which the particle density becomes 
less than 1/fm 3 . 
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FIG. 1. The logarithm L (see section V A) as a function of time (in fm) in the saturation scenario 
at RHIC and LHC. 
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FIG. 2. Various observables as functions of time in the saturation scenario. The upper panel 
shows the time evolution of the energy density and of the longitudinal and transverse pressures 
(times 3) at RHIC. The second and third panels show the evolution of the average longitudinal and 
transverse squared momenta per particle (pi) and (p^_) at RHIC and LHC respectively. Comparing 
with the curves obtained in [8], one observes a semi-quantitative agreement (in the above figure, 
everything is expressed in units of GeV, as in [8]). 
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FIG. 3. Time evolution of the ratios R^, defined in Eq. ( p. 4 ) in the saturation scenario at 
RHIC (left) and LHC (right). The quantities 1 — Rk give a measure of the anisotropy of the 
microscopic distribution f(p, i) for different momentum scales, the larger k the lower the probed 
momentum scale . These curves are obtained with as = 0.3 and the logarithm L = ln((pj) / m^) 



(see section VA) 
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FIG. 4. Time evolution of the ratios Rk, defined in Eq. ( p. 4 ) in the minijet scenario at RHIC (left) 
and LHC (r ight) . These curves are obtained with as = 0.3 and the logarithm L — ln((pf)/7n|>) 
(see section VA). 
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FIG. 5. Sensitivity of the results (here for the time evolution of the ratio Ri = Pl/Pt) with 
respect to the details of the description in the saturation (upper panels) and minijet (lower panels) 
scenarios at RHIC (left panels) and LHC (right panels) energies. Compared to the results shown 
in Figs. and ^, the dotted curves are obtained by replacing the logarithm L — » 2L and the dashed 
curves by replacing L — > L/2. The solid curves in th e minijet scenario correspond to the case where 
the coupling constant runs with time (see Eq. (5.4)). 
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